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ABSTRACT 

Recent measurements of the Zeeman effect in dark-cloud cores provide important tests for theories of cloud 
dynamics and prestellar core formation. In this Letter we report results of simulated Zeeman measurements, 
based on radiative transfer calculations through a snapshot of a simulation of supersonic and super-Alfvenic 
turbulence. We have previously shown that the same simulation yields a relative mass-to-flux ratio (core versus 
envelope) in agreement with the observations (and in contradiction with the ambipolar-drift model of core 
formation). Here we show that the mass-to-flux and turbulent-to-magnetic-energy ratios in the simulated cores 
agree with observed values as well. The mean magnetic field strength in the simulation is very low, B = 0.34 jiG, 
presumably lower than the mean field in molecular clouds. Nonetheless, high magnetic field values are found in 
dense cores, in agreement with the observations (the rms field, amplified by the turbulence, is B lms = 3.05 /iG). 
We conclude that a strong large-scale mean magnetic field is not required by Zeeman effect measurements to 
date, although it is not ruled out by this work. 

Subject headings: ISM: magnetic fields — stars: formation — MHD — radiative transfer 



1. INTRODUCTION 

Supersonic turbulence can naturally lead to the complex 
structure and the large density contrast found in molecular 
clouds, a process often referred to as turbulent fragmenta- 
tion. Magnetic fields can affect the outcome of the turbulent 
fragmentation by reducing the density contrast of shocks, by 
providing support against gravitational collapse, and by en- 
hancing angular momentum transfer. To test and calibrate 
star-formation models and simulations based on observational 
data we must determine the magnetic field strength in the ob- 
served star-forming regions. This can be achieved by mea- 
surements of the Zeeman effect on molecular emission lines, 
but only few such measurements have been obtained in dark 
cloud cores to date (32 measurements and 9 detections in 
iTroland & Crutcherl (120081) ). 

The detection of the Zeeman effect on a handful of molec- 
ular cloud cores is not sufficient to provide a direct estimate 
of the mean or mean-squared magnetic field in star-forming 
regions, as that would require measuring the field strength at 
many positions in the same molecular cloud, including low 
density regions. Furthermore, as discussed in §3, our radiative 
transfer calculations show that the OH Zeeman effect gives a 
magnetic field corresponding to the average value of p a B\ os 
along the line of sight, with a « 1.5. This is because the 
OH emission in the 1665 MHz and 1667 MHz lines that are 
used for measuring the magnetic field strength increases with 
both gas density and excitation temperature, and the excita- 
tion temperatures increase with density until the transitions 
become thermalized at n sa 10 3 cm" 3 . Therefore, even in the 
rare locations where the Zeeman effect is detected, it does not 
provide a direct estimate of the mean field in that line of sight, 
but only the field strength in the densest gas. However, it may 
be possible to constrain the mean magnetic field strength us- 
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ing numerical simulations, as long as observable quantities 
derived from the simulations are found to depend on the mean 

m agnetic field. 

Pad oan & Nordluridl(ll999l) discussed several observational 
tests of numerical simulations of supersonic MHD turbulence 
that could be used to constrain the mean magnetic field. They 
concluded that the mean magnetic field in molecular clouds 
is significantly weaker than previously assumed in models of 
molecular clouds and theories of star formation, suggesting 
that molecular-cloud turbulence is super-Alfvenic on the aver- 
age, meaning on scales of few to several parsecs. This super- 
Alfvenic model of star-forming regions was recently used to 
generate simulate d measurements of the Zeeman effect on 
18 cm OH lines (iLunttila et al.ll2008l) . It was shown that a 
super-Alfvenic turbulence simulation with the characteristic 
size, density, and velocity dispersion of star-forming regions 
could produce dense cores w ith the same |ffio S |-jV relation as 
observed cores. Furthermore. ILunttila et al.l (120081) computed 
the relative mass-to-flux ratio IZ^ , defined as the mass-to-flux 
ratio of the core divided by that of the envelope, following th e 
observational procedure proposed by Crut cher et al.l (2009). 
They found a large scatter in the value of IZ^, and an aver- 
age value of IZft < 1 , in contrast to the ambipolar-drift model 
of core formation, where the mean magnetic field is stronger 
and only 7£ u > 1 is allowed. The observational results of 
ICrutcher et al.l (120091) confirmed TZq < 1 in observed cores, 
as predicted by ILunttila et all d2008l) for the super-Alfvenic 
model. 

In this Letter we present further evidence that the same 
super-Alfvenic simulation compares well with the observa- 
tional data. We use simulated OH Zeeman measurements to 
compute the mass-to-flux ratio relative to the critical one, A, 
and the ratio of turbulent to magnetic energies, /3 tU rb, in molec- 
ular cores selected from simulat ed maps. We follow closel y 
the observational procedure of ITroland & Crutcher! (2008), 
and find mean values of A and /3 tm -b in good agreement with 
their observational results. 
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2. NUMERICAL SIMULATION OF SUPER-ALFVENIC TURBULENCE 

This work is based on the same simulation of supersonic 
a nd super-Alfvenic tu rbulence (and the same snapshot) used 
in lLunttila et al.l (120081) . The simula tion was run on a m esh of 
1000 3 zones with the Stagger Code dPadoan et al.ll2007l) . with 
periodic boundary conditions, isothermal equation of state, 
random forcing in Fourier space at wavenumbers 1 < k < 2 
(k = 1 corresponds to the computational box size), uniform 
initial density and magnetic field, and random initial velocity 
field with power only at wavenumbers 1 < k < 2. The rms 
sonic Mach number is Ai s = ct v .3d/c s = 8.91. 

The initial value of the ratio of gas to magnetic pressure 
in the simulation is = {pc\)/(Etf j /'&ir) = 22.2, where cs is 
the sound speed, and Bo is the initial (uniform) magnetic field 
strength. At the time corresponding to the snapshot used in 
this work, the rms magnetic field strength has been amplified 
by the turbulence, giving a value of j3 = 0.2, defined with the 
rms magnetic pressure. This corresponds to an rms Alfvenic 
Mach number of M- d = (/3/2) 1/2 cr v j D /c s = 2.8, so the turbu- 
lence is super-Alfvenic also with respect to the rms Alfven 
velocity. With respect to the Alfven velocity corresponding to 
the mean magnetic field (/3; = 22.2), the rms Alfvenic Mach 
number is much larger, Ai- d .i = (/3i/2) 1 / 2 tr V 3D/c s = 29.7. Val- 
ues of parameters scaled to physical units are given in the next 
section, as the radiative transfer calculations require physical 
values of size, temperature, and mean density. 

3. SIMULATED ZEEMAN EFFECT 

For the computation of synthetic Zeeman spectra the data 
cube is scaled to physical units. The size of the grid is fixed 
to L = 9 pc, the mean densit y to (njfy)) = 67 cm" 3 ( typical for 
that scale in the sample of Falgar one et al.l dl992h ). and the 
kinetic temperature to Jkin_= 10 K. With this scaling, the mean 
magnetic field is (B Xl B y ,B z ) = (0.0,0.0,0.34) /iG and the rms 
field Z? rms = 3.05 /xG (much larger than the mean magnetic 
field, due to the turbulent amplification mentioned above). We 
assume a constant fractional OH abundance of [OH]/[H] = 
4.0 x 10" 8 from lCrutcheri ( [19791) . 

We simulate Zeeman splitting observations of 1665 and 
1667 MHz OH lines that are commonly used for measuring 
magnetic field strengths in molecular clouds. Full radiative 
transfer cal culations with our line radiative transfer program 
(lJuvelall 19971) . using a datacube resampled to a resolution of 
256 3 cells, showed that the cloud is optically thin and the ra- 
diation field is approximately constant throughout the cloud. 
Thus, the level populations only depend on the local density. 
Comparison with our full radiative transfer calculations indi- 
cate that the errors in the excitation temperatures of the 1665 
and 1667 MHz OH transitions are < 1 K. Using the estimated 
level populations, the coupled radiative transfer equations for 
the four Stokes parameters are integrated along the line of 
sight. We use a resolution corresponding to 512 3 computa- 
tional cells for computing the synthetic Zeeman spectra. To 
increase the number of simulated observations, the calcula- 
tions are carried out for the x, y, and z directions. 

To compare the results with observations, we simulate ob- 
servations at cloud distances of D = 150 pc, D = 300 pc, and 
D = 1000 pc, yielding total angular sizes of ~ 3.4°, ~ 1.7°, 
and ~ 31'. We compute position-position-velocity data cubes 
by simulating observations of the 1665 MHz line. The syn- 
thetic observations are made with a 3' (fwhm) beam, corre- 
sponding to the angular resolution of the Arecibo telescope, 
and using a channel separation of 0.05 km s" 1 . The left panel 



of Figure 1 shows the simulated 1665 MHz OH integrated in- 
tensity map (for the y direction), assuming a cloud distance of 
D = 150 pc. The corresponding OH peak brightness tempera- 
ture map is shown in the right panel of Figure 1 . 

Figure 2 compares the value of Blos estimated from the 
Zeeman effect (right panel) with the actual mean value of the 
line-of-sight magnetic field computed directly from the three- 
dimensional data cube, (BLOS.tme) (left panel), for the same 
map shown in Figure 1 . The Zeeman effect tends to select the 
densest regions along the line of sight, where the magnetic 
field is stronger (dense cores originate from super-Alfvenic 
compressions that amplify the magnetic field components on 
the plane perpendicular to the direction of compression). As a 
result, the Zeeman effect greatly overestimates the mean mag- 
netic field on lines of sight passing through dense cores - in 
the y-direction map, |(-SLOS,true)| reaches a maximum value 
of only 3.7 /iG, while |#los| reaches a maximum value of 
48 /iG. Because the Zeeman effect is weighted by OH emis- 
sion, hence by both density and excitation temperature (which 
increases with density up to n sa 1000 cm" 3 ), the estimated 
magnetic field strength corresponds approximately to the line- 
of-sight mean value of //'.Blos, with a ~ 1.5. The difference 
between the line-of-sight mean magnetic field and its value es- 
timated through the Zeeman effect is further illustrated in the 
left panel of Figure 3, showing a scatter plot of Blos (from 
the Zeeman effect) versus (BLOS.true)- The right panel of Fig- 
ure 3 shows that Z?los is well correlated with (p L5 BLOS,trae), 
hence the Zeeman effect is strongly biased towards the mag- 
netic field strength of the densest regions along the line of 
sight. 

We use the computed OH emission position-position- 
velocity dat a cubes to select dens e cores with the clumpfind 
algorithm dWilliams et alj 1 1995b . Before applying the 
clumpfind routine, the data cubes are resampled to an angu- 
lar resolution of ^ 1 .2' (approximately Nyquist sampled), 
and uncorrelated Gaussian noise with rms of 0.08 K is added 
to simulate observational noise. The clumpfind threshold and 
stepsize parameters are both set to 0.4 K (5cr rms ). The algo- 
rithm provides a list of detected cores a nd their sizes , which 
are needed in the subsequent analysis. iPineda et al.l (2009) 
have recently shown that the clumpfind algorithm is sensi- 
tive to the values of the parameters. Running the clumpfind 
with other parameter values would yield a set of cores with 
different properties. Nevertheless, our choice of parameters 
produces a repres entative ensemble of cores similar to the ob- 
served sample of iTroland & Crutch er (2008). The core posi- 
tions are shown as circles on the maps of Figure 1 , with the 
circle size equal to the size of the 3' (fwhm) beam. The cores 
selected with the clumpfind algorithm do not match very well 
the column density structure revealed by the integrated inten- 
sity map. This is to be expected because the clumpfind algo- 
rithm tends to isolate real cores, rather than structures in pro- 
jection, by selecting distinct features in radial velocity space. 
As shown by the right panel of Figure 1, all the cores corre- 
sponds to local maxima in the OH peak brightness tempera- 
ture map. 

Each of the detected cores is observed with a 3' beam, tar- 
geting the position of peak intensity in the core. The sim- 
ulated observations consist of Stokes I and V spectra of the 
1665 and 1667 MHz lines for each core. To simulate noise 
in the observations, uncorrelated Gaussian noise is added to 
both I and V spectra. The line-of-sight magnetic field is deter- 
mined by least-squares fitting the numerical derivative of the 
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FIG. 1 . — Left: Simulated 1665 MHz OH integrated intensity map in the y direction, convolved with a 3' beam, assuming the model cloud is at a distance 
of 150 pc. Circles show the size of the beam (fwhm) and the positions of the cores. Right: Core positions shown on the corresponding OH peak brightness 
temperature map. 
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FIG. 2. — Left: Mean value of the line-of-sight magnetic field computed directly from the three-dimensional data cube. Right: Line-of-sight magnetic field 
estimated from the simulated Zeeman measurements. 



Stokes I spectrum to the Stokes V spectrum, as it is usually 
done with Zeeman splitti ng observations (e.g. ICrutcher et all 
ll993l;lBourke et ail uOOl). The magnetic field strength is de- 
termined separately from both the 1665 MHz and the 1667 
MHz lines, and the analysis uses the average of the two mea- 
surements weighted by their inverse squared errors in B. 

Column densities are estimated from the simulated spec- 
tra. Assuming that the line is optically thin, the column den- 
sity of OH [cm" 2 ] is obtained as Af(OH) = a(l - T^/T^W, 
where W is the integrated line area [Kkm s -1 ], 7t, g is the back- 
ground continuum brightness temperature (2.73 K), T ex is the 
excitation temperature of the transition, and a is 4.04 x 10 14 
for t he 1665 MHz l ine and 2.24 x 10 14 for the 1667 MHz 
line dCrutcherl [19791) . We use T ex = 10 K in the analysis 
of our synthetic observations. The total H2 column density 
is obtained as N(H 2 ) = jV(OH)/8 x 10" 8 . Core masses and 
vo lume densities are calcula ted with the same formulas as 
in iTroland & Crutcherl d2008l) : n(H 2 ) = JV(H 2 )/2r, M(OH) = 
7rr 2 Af(H 2 )2.8/«H, where r is the core radius and toh is the mass 
of an H atom. 

4. RESULTS 

In the synthetic maps, the clumpfind algorithm finds a to- 
tal of 105 cores at D = 150 pc, 40 at D = 300 pc and four at 
D = 1000 pc. Our simulated observations of the cores use a 
noise level corresponding to a 30-hour integration with a sys- 



tem temperature of 3 5 K, similar to parameters in recent OH 
Zeem an surveys (e.g. lBourke et al .1200 U ITroland & Crutcherl 
2008). As in actual Zeeman splitting observations, due to 
very low amplitude of the Stokes V spectrum, statistically sig- 
nificant detections of the magnetic field are rare. With the 
commonly used criterion, |Blos| > 3o£(los), there are only 
24 detections at D = 150, four at D = 300 pc and no detec- 
tions at 1000 pc. The ratio of the number of 3a detections to 
the number of observed cores, 28/149, i s similar to what has 
been achieved in Zeeman surveys (e.g. JCrutcher et af1ll993t: 
iBourke et alj200U ITroland & Crutcherll2008l) . 

In ITroland & Crutc her (2008), the authors employ an addi- 
tional criterion for determining whether the measurement is 
a detection: 1665 MHz and 1667 MHz lines must yield the 
same Blos within the measurement uncerta i nties. If the same 
threshold is used as in ITroland & Crutc her (2008), AB L0S = 
|fl L os(1665)-fi L qs(1667 )J < 1.9a A B, three of the 2 8 3cr de- 
tections fail this criterion. ITroland & Crutchen (120081) also use 
a subjective criterion that the V spectrum must look consis- 
tent with detection of the Zeeman effect. This is the case for 
all 28 3cr detections in our simulations. However, in some 
of th ese spectra (and in some of the observed spec tra as well 
(e.g. lBourke et al.ll200U ITroland & Crutcherll2008l) ) there are 
signs of multiple velocity components with differing line-of- 
sight magnetic fields. 
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FIG. 3. — Left: Isodensity contours of the line-of-sight magnetic field estimated from the Zeeman measurements on the 1665 MHz OH line, versus the mean of 
the true line-of-sight magnetic field. Right: Isodensity contours of the line-of-sight magnetic field estimated from the Zeeman measurements versus the average 
of the true line-of-sight magnetic field multiplied by p 5 . The contour levels in both panels are 40, 400, 4000, and 40000 pixels fiG~ 2 . 



TABLE 1 

Physical parameters from the synthetic observations. 



D < 200 pc 


200 pc< D < 400 pc 


D > 400 pc 


simulations T&C08 1 


simulations T&C08 


simulations T&C08 


n=W5 n=l6 


n = 40 ;i = 13 


n=4 n=5 



|B L os| a 6.89±0.65 6.54±1.39 7.35±1.15 11.1 ±2.7 11.2±4.7 17.2± 10.8 

r b 0.206±0.010 0.234±0.015 0.413±0.025 0.431±0.042 0.751±0.054 1.53±0.22 

AV C 0.498±0.021 0.622±0.063 0.856±0.072 0.890±0.110 1.22±0.20 1.67±0.10 

W(H 2 ) d 3.73±0.16 4.32±0.56 3.61 ±0.21 3.78±0.35 3.29±0.44 6.58±1.00 

n(H 2 ) e 3.22±0.13 3.09±0.40 1.53±0.10 1.48±0.13 0.707±0.072 0.765±0.125 

M(OH) f 16.5±2.2 19.4±4.6 57.1 ±9.1 59.0±12.0 139±33 1090±210 

M(virial) f 16.9 ±2.2 28.0 ±7.2 105 ±21 118 ±30 316 ±108 1070 ±160 

A,A TC s 4.12|| 5.04|| 3.76«4 2 . 63 4.72 2 . 30 8jt 3 . 30 15.5 

Art, A U rb.TC g 1.60?$ 2.78?;* 2.14?$ 1.07«J 0.9l£g 0.97*1° 



1 Observations from Troland & Crutcher 1 2008). a Mean line-of-sight magnetic field, /iG. b Mean core radius, pc. c Mean fwhm, km s '. d Mean column density, 
10 21 cm -2 . 1 Mean volume density, 10 3 cm -3 .' Mean mass, Mq. 8 median 025quantne 



iTroland & Crutcher] d2008l) observed cores at distances 
ranging from 140 pc to 2000 pc. Because the 3' beam se- 
lects structures at very different scales at the extremes of 
this distance range, we compare the simulation results at 
each assumed distance only to the observations of cores 

at similar distance s. Thus, we consider separately the 

ITroland & Crutcher] d2008l) observations of cores at D < 200 
pc, 200< D < 400 pc, and D > 400 pc. The mean values of 
physical parameters of the cores as determined from the syn- 
thetic observations are given in Table 1 . For comparison, we 
also list the corresponding values calculated from the cores in 
ITroland & Crutcher] d2008l) . The results for key physical pa- 
rameters determined from ou r synthetic observations are con- 
sistent with the results from Trol and & Crutcher] d2008l) for 
th e two closest distances. Th e most distant cores observed 
bv lTroland & Crutcher! d2008) have radii approximately twice 
larger than the cores in our simulations; they are also approxi- 
mately 8 times more massive and their linewidths are approx- 
imately \pl ti mes larger than the s i mulate d cores. 

Following ITroland & Crutchen d2008l) . we define A = 
(M/$) b s /(M/$) cl -i t , the ratio of observed and critical mass- 
to-flux ratios (above the critical mass-to-flux ratio the mag- 
netic field cannot support the cloud against gravitational 
collapse). Th e theoretically determined critical mass-to- 
flux ratio from lNakano & Nakamural (1 19781) yields A = 7.6 x 
10" 21 /Y(H 2 )/|fi L os|, where N(U 2 ) is in cm -2 and B h0S is in 



fiG. Because only the line-of-sight component of the mag- 
netic field is known, values of A are not meaningful for in- 
dividual clouds. Instead, the values given in Table 1 are cal- 
culated based on the mean values for /V(H2) and |Blos|- Be- 
cause A is a nonlinear function of the observed parameters, its 
value and uncertainty are estimated using Monte Carlo sam- 
pling. The resulting error distributi on has a long tail towards 
high values of A. Furthermore, Troland & Crut cher] (1200 8l) 
define the ratio of turbulent and magnetic energies, /3 tU rb = 
Eturb/Emag, as /3 t ui b = 0.32»(H2)AV^ T /( 3g LOS ), where »(H 2 ) 
is in cm" 3 , fi LO s is in ^G, and AV NT = VAV 2 - 0.027 km 2 s" 2 
is the non-thermal linewidth in km s" 1 . As with the case of A, 
instead of calculating the ratio for individual clouds, we use 
the mean values of n(H2), AV, and Blos, and estimate the er- 
ror distribution with Monte Carlo sampling. The values of A, 
Atc> Aurb, and /? tur b,TC in Table 1 are the median values com- 
puted from the Monte Carlo sampling around the mean values 
of the observed parameters. 

The mean values of all the observed and derived parameters 
from t he simulation are consisten t with those from the sam- 
ple of ITroland & Crutcher] (120081) . As a result, the values of 
A from the simulation, A, are also consi stent with the corre- 
spond ing values from the observations of lTroland & Crutcher! 
(2008), Axe- The value of A is slightly smaller than Atc for 
cores at distances D < 200 pc and D > 400 pc, and slightly 
larger than Atc for cores at distances in the range 200 pc 
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D=150 pc 




FIG. 4. — - Left: Inferred line-of-sigh t magnetic field as a function of H2 column density. Black symbols show the results from our simulation, red symbols are 
cores from Troland & Crutcher 1 2008). Squares show 3<r detections with ±lcr errorbars, and triangles are no n-detections. The black solid line shows the mean 
mass-to-fiux ratio, A, and the red solid line the mean mass-to-flux ratio calculated from Troland & Crutcher (2008) observations, Atc- The dashed line is the 
critical mass-to-fiux ratio with no geometrical correction, A = 1 . Right: Magnetic energy versus turbulent kinetic energy for the same cores and with the same 
symbols as in the left panel. The dashed line corresponds to /3 lur (, = Eturb/finag = 1- The solid black line shows the mean energy ratio, /3 tU rb> an d the red solid 
line the mean energy ratio calculated from Troland & Crutcher (2008) observations, Aurb.TC- The upper panels compare the results of simulated observations at 
D = 150 pc with observations of cores at D < 200 pc. The lower panels are for simulated cores at D = 300 pc and observations in the range 200 pc < D < 400 pc. 

< D < 400 pc. The same applies to the comparison of the 
mean energy ratios_from the simulation, /3 tur b, with that from 
the observations, /3 tU rb,TC- In the case of (3 tur b, however, the 
differences are larger because this ratio depends on the square 
of the measured linewidth and magnetic field strength. The 
values of |Blos| and iV(H2) are shown in the left panels of 
Figure 4 for all cores at distances D < 200 pc (upper panel) 
and 200 pc < D < 400 pc (lower panel). The right panels 
of that figure show the values of E mag = 47rr 3 Z?Los/(^ 7r ) an d 
£ turb = (3 /2)M(OH)[AVnt/(2 ^2 log 2)] 2 for the same cores. 



In both plots, the cores selected from the simulation (black 
symbols) cover approximately the same regions of parameter 
space as the observed cores (red symbols). 

5. CONCLUSIONS 

This work shows that supersonic and super-Alfvenic turbu- 
lence generates dense cores with physical properties similar 
to those of observed cores. The mean values of the mass-to- 
flux ratio and of the turbulent-to-magnetic-energy ratio mea- 
sured in the simulated cores are c onsistent with the obse rved 
values. As previously shown by ICrutcher et alj d2009t) . the 
observed relative mass-to-flux ratio (cores versus envelope) 
is consistent with the prediction of super- Alfvenic turbulence 
presented in Luntti laet al.l d2008l) as well. The simulation 
used in this work has a very low mean magnetic field strength, 
B = 0.34 /j,G, presumably lower than in real molecular clouds. 
We have chosen this extreme example of super-Alfvenic tur- 
bulence to show that the observed core magnetic properties 



do not require a large mean magnetic field on the large scale. 
The reason why a field strength much larger than the mean 
may be found in dense cores, is that the field is locally ampli- 
fied by the turbulence, particularly by the strong compressive 
motions creating the cores. 

We conclude that all Zeeman measurements to date are con- 
sistent with the super-Alfvenic model of molecular cloud tur- 
bulenc e and core for mation proposed bvlPadoan & Nordlundl 
d!999t) . As shown bv lCrutcher et alj (12009). the observations 
also appear to contradict a model of core formation based on 
large-scale magnetic support and quasi-static evolution con- 
trolled by ambipolar-drift. While we have proved that a strong 
large-scale mean magnetic field is not required by the Zeeman 
effect measurements, a strong mean field cannot yet be ruled 
out in the context of supersonic turbulence. In the presence 
of supersonic turbulence, however, ambipolar drift is not the 
dominant mechanism of core formation, because cores can be 
formed dynamically by random compressions along the field. 
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